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METHOD OF OPTIMIZING PARAMETER VALUES IN A PROCESS OF PRODUCING 

A PRODUCT 

BACKGROUND OF THE INVENTION 

Field of th invention 

5 The invention relates to the process optimization field, and more particularly to a 

method of optimizing parameter values in a process of producing a product which is 
characterized by properties affected by the selected parameter values. This invention is 
applicable in different industries, such as the pharmaceutical, chemical, cosmetics, 
plastics, petrochemical, agriculture, metallurgy and food industries, as well as many 
10 other commercial and industrial applications. 

Description of prior art 
Processes for production of complex compositions such as those found in many 
m pharmaceutical products generally require the mixing of many ingredients according to 
rf specific process parameters regarding formulation and production technology, to provide 
\Fh5 the product with properties at a level offering satisfactory performance according to 
If predetermined specifications. In such complex production processes, it is not unusual 
y that some process parameters involved exhibit interfering effects on the desired 
« properties, further complicating the process design. Where possible, the designer may 
[q try to adapt the set of process parameters from known data derived from previous similar 
1^20 processes, and/or rely on conventional trial-and-error experimental schemes to optimize 
0 the set of process parameters values, in order to meet the product specifications. 
• However, as the processes become more complex, optimization in such 
multidimensional space with high accuracy requirements turns out to be an extremely 
difficult task, even for the highly skilled designer. That limitation is particularly 
25 problematic in the design of pharmaceutical products, where one or more active 
substances mixed with a variety of excipients (e.g. carriers) must be produced in the 
form of a stable and highly effective standard delivery system such as a tablet, capsule, 
suspension, cream or injection, or even controlled release systems such as skin carriers 
and implants. 

30 In the past years, many techniques have been developed to assist the process 

designer or formulator in optimizing values of parameters governing processes. These 
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techniques aim at quantify existing relations between parameters and associated 
desired product performance characteristics. A conventional technique known as the 
Full Factorial Matrix (FFM) method consists of statistically deriving a behavior relations 
for the properties from a set of experimental runs of the process using selected initial 
5 values for the parameters. The established model being generally nonlinear, optimized 
parameter values are then derived using an optimization method such as the 
Multisimplex method described in "Practical Methods of Optimization" J, Wiley & Sons, 
Chichester, 2d, (1987), which essentially consists of linearizing the behavior functions 
related to the parameters according to straight lines or planes of different random 
10 directions. For any given property behavior relation of n parameters to be optimized in 
order to either minimize or maximize that behavior relation with or without constraints on 
the parameter values, a recursive estimation of the property is then performed using an 
m initial set of parameter values according to a selected direction, until the obtained value 
for the property does not significantly vary in that direction. Then, a last unfavorable set 
Ul5 of parameters is used as a new starting point for a following recursive estimation 
according to a different direction. Successive recursive estimation steps are performed 
;rj until the resulting value for the property no longer significantly vary in any new direction. 
« When applied to a model comprising a plurality of property behavior relations, the 
Iq Multisimplex method allows a unique objective function to be created by proper 
1^0 transformation of the relations to adapt to different scales and/or units and by 
ifl ? associating a relative importance weight to each property, either subjectively or through 
'** • fuzzy logic algorithms. 

The known optimization processes based on Full Factorial Matrix-Multisimplex 
methods suffer from several drawbacks. As a general rule, the number of experimental 
25 runs required to obtain a model of sufficient reliability is proportional to the total number 
of significant parameters involved. Therefore, the cost and time frame of the 
experimental work will therefore be essentially proportional to the number of runs 
required. Although a variant of the method known as the Fractional Factorial Matrix has 
been proposed in order to reduce the number of runs to be performed, the provided 
30 reduction of experimental runs may not significantly reduce the total cost and time frame 
of the work required to complete the design of a complex product involving many 
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production technologies. While adequate formulations complying with constraints 
imposed on the parameter values can nevertheless be obtained, these formulations 
generally cannot be qualified as optimal when comparing actual property performance 
with desired property values set forth in the product specifications. 
5 A technique which attempts to improve parameter optimization in process design 

is disclosed in European Patent Office laid-open patent application publication number 
0,430,753 dated June 5, 1991 and in US patent No. 5,218,526 issued on June 8, 1993 
to Mozzo. According to the technique in Mozzo, from a set of property relations 
expressed in terms of parameters which is obtained by standard statistical methods 
10 using the results of a number of experimental runs of the process, a corresponding set 
of property relations expressed in terms of weighted parameters is derived. For each 
actual value of a parameter, a first weighting is expressed as the ratio of: (a) the 
gi deviation of the actual value from the mean value of the parameter over the 
experimental range, on (b) the range between extreme values for that parameter over 

Uli5 the experimental range. Then, a goal function is established in term of deviations 

IV 

,p between weighted values of property values as estimated by the property relations and 
corresponding weighted values of specified goal values for the properties. For each goal 
« value of a property, a second weighting is expressed as the ratio of: (a) the deviation of 
,g the actual value from the mean value of the property over the experimental range, on (b) 
!~20^ the range between extreme values for that property over the experimental range. Then, 
\Q ' according to a recursive geometric algorithm aimed at successively minimizing the 
established goal function, a set of optimal parameter values is generated. While being 
an improvement over the conventional Full/Fractional Factorial Matrix - Multisimplex 
methods regarding the capability to consider specified goal values for the properties, the 
25 weightings as taught by Mozzo do not reflect the relative importance of the properties 
involved, and that limitation may therefore affect the convergence of the algorithm 
toward an optimal solution. 

A review of modem techniques and software systems for the design of 
pharmaceutical product formulations is given in "Intelligent Software System For 
30 Pharmaceutical Product Formulation" R.C. Rowe, Pharmaceutical Technology, March 
1997. In that paper, expert systems, rule induction algorithms, case-based reasoning 
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algorithms, neural networks and genetic networks are presented as modern tools for 
supporting formulation design, and a number of available software systems using some 
of these tools are summarized. As indicated in the Rowe paper, although a 
knowledgeable expert system could be a powerful tool to assist the process designer in 
5 the formulation task, its development is generally a high risk, time consuming and 
expensive process. Rule induction is a knowledge-based algorithm which allows 
hierarchical classification of objects, using statistical methods which are found generally 
effective only if the input data is continuous, which is often not the case in practice. 
Moreover, since rule induction is limited to establishing whether or not a given object is 
10 close to another, it generally cannot provide an optimal solution. Case-based reasoning 
is a knowledge-based iterative technique which can be used to design formulations, 
which consists of matching the desired specifications for the product with the 
m specifications of the most relevant known formulation(s), and adapting the selected 
)^ formulation(s) as necessary, followed by an evaluation. Although effective for optimizing 
Ml 5 the parameters of a variant process from a family of similar processes and 
,p corresponding formulations, case-based reasoning generally cannot be used where the 
;y design of a significantly different formulation is contemplated. As to neural networks, in 
which each neuron input is modified by a weight associated with that neuron, they 
2 appear to be effective tools for assisting formulation design only in cases where no 
:=g0 constraint applies on either the parameter or property values, such cases being rarely 
\Q found in practice. Finally, regarding the genetics algorithms, they are cyclic methods 
~ based on Markov chains for predicting from a starting point a solution likely to result from 
a sequence of operations, in order to allow making changes to obtain a desired solution. 
Since these changes are generally made arbitrarily, in most cases, the resulting solution 
25 cannot be considered as optimal. 

SUMMARY OF THE INVENTION 
It is therefore an object of the present invention to provide a systematic method of 
optimizing parameter values in a process for producing a product which minimizes the 
number of experimental runs required to obtain an optimal solution complying with the 
30 product specifications. 
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According to the above object, from a broad aspect of the present invention, there 
is provided a method of optimizing parameter values in a process of producing a 
product, the process being essentially controlled by a set of n parameters X { affecting 
a set of k properties Y s characterizing the product. The method comprises the steps of: 
5 i) assigning values to a set of k property weights w y representing relative importance of 
the properties Y } for the characterization of the product; ii) establishing property 
behavior mathematical relations giving an estimated property Ye } for each property 7. in 
terms of the parameters X { from given parameter data and associated property data; iii) 
using the property weights w y to establish a goal function in terms of property weighted 
10 deviations between the estimated properties Ye s and corresponding specified goal 
m values for the properties Y J ; and iv) optimizing the goal function to generate a set of n 

a 

$•* optimal parameter values for the parameters X. . 
Ul 

ry According to a further broad aspect of the present invention, there is provided a 

jn method of producing a pharmaceutical product using optimized process parameter 
"S5 values, the process being essentially controlled by a set of n parameters X s 
Q characterizing a formulation for the product, the parameters X £ affecting a set of k 
properties Yj characterizing the product. The method comprises the steps of: a) 
conducting a number of / of experimental runs of the process each using a selected 
distinct set of values for the parameters X i covering substantially all extreme values 
20 within a chosen range of values for each one of the parameters X i , wherein / is at least 
equal to n + 1 and is substantially less than a number used in the Fractional Factorial 
Matrix method; b) measuring values for the properties Y } characterizing the product in 
each of the / experimental runs, whereby parameter data and associated property data 
are obtained from the selected distinct set of values for the parameters X i and the 
25 measured values for the properties Y s , respectively; c) determining an importance of the 
properties Y- for the characterization of the product, comparing the importance of the 
properties Yj relative to one another, and assigning values to a set of k property 
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weights w. representing a relative importance of the properties Y. for the 
characterization of the product; d)calculating a set of optima! parameter values for the 
parameters X t using the measured values for the properties Y j and the assigned values 

of the set of k property weights Wj ; and e) producing the pharmaceutical product using 
5 the optimized process parameter values X t calculated in the previous step. 

BRIEF DESCRIPTION OF THE DRAWING 
The invention will be better understood by way of the following detailed 
description of a preferred embodiment with reference to the appended drawings, in 
which: 

10 Fig. 1 is a block diagram of a software system that can be used to carry out the 

method according to the present invention according to the preferred embodiment; and 
(II Fig. 2 is a flow chart representing the preferred embodiment of the method 

|J according to the present invention. 

U1 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

J3 5 In the following description, a preferred embodiment of the present invention 

Ig applied to product formulation design will be described. However, it is to be understood 
* that the present invention can be also be used to optimize parameter values of 
J© processes related to the production of many types of products which cannot be 
R associated with a formulation, while being characterized by a number of properties 

^0 affected by process parameters, such as biotechnological products, electronic 

•0 

~" components, etc. 

Referring now to Fig. 1, there is generally designated at 10 a computer system 
which is programmed to carry out a method according to the present invention. The 
system preferably comprises a knowledge base 12 where prior formulation/process data 

25 and competing products data are stored. For the purpose of pharmaceutical formulation 
design, knowledge base 12 contains process data related to ingredient proportions, 
experimental conditions and results over time, production technologies used, etc. The 
system 10 further comprises a property weighting module 14 which generates a weight 
value for each one of a number k of identified properties according to an initial modeling 

30 of the problem and property comparison data presented to the module 14. System 10 



SOR 14081-1 USPR 



6 



further comprises an evaluation module 18 fed by the property weights generated by 
module 14, to generate a global relative importance vector of dimension [k] for the k 
properties. System 10 is provided with an experimental data entry module 16 through 
which property values obtained from experimental runs using different sets of parameter 
5 values for the process can be entered and stored for later use by several modules of 
system 10. Linked to receive data from modules 12, 14 and 16, is an evaluation module 
18 which can generates a ranking of the sets of property values selected from the 
knowledge base and the optimal set of parameter values obtained through optimization. 
System 10 further comprises a parameters reduction module 22 to retain only those 
10 parameters having a significant effect on the considered properties. Module 22 is 
particularly useful in cases where the number of parameters involved is relatively large, 
usually greater than 8 where a computer provided with a standard high-performance 
«j microprocessor is used. The S-Plus™ statistical software from MathSoft may be used in 
U module 22 to carry out the Stepwise method to select the variables . System 10 is 
pjs further provided with a parameter interaction module 20, the function of which consists 
*P of identifying by statistical analysis from experimental data, which remaining parameters 
si are significantly correlated. The S-PIus statistical software from MathSoft can also be 
IU used to program module 20 in which the appropriated correlation methods are applied to 

*S the data. It is to be understood that module 20 is unnecessary where all parameters are 

Q 

\g0 independent one another. 

0 Modules 16, 20 and 22 are linked to a property behavior models module 24 that 

uses experimental data, parameter interaction data and remaining significant 
parameters for determining an optimal mathematical model for each property which is 
likely to better estimate that property. The model data as generated at module 24 is fed 
25 to a property behavior relation module 26 that also receives experimental data from 
module 16 to statistically estimate polynomial coefficients to be incorporated within the 
established property behavior models, thereby generating a behavior relation for each 
property. The S-Plus statistical software from MathSoft may be used to program module 
26 to apply the appropriate regression methods to the data.System 10 is further 
30 provided with a goal function module 28 linked to property weighting module 14 and 
property behavior relation module 26 to generate, from specified goal values for the 
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properties, a goal function in term of property weighted deviations between properties as 
estimated by the behavior relations and the corresponding specified goal values for 
these properties. 

An optimization module 30 is provided to optimize the goal function as 
established by module 28 through successive iterations and according to the type of 
each variable (discrete or continuous) and according to one or more ranges specified as 
constraints imposed on one or more optimal parameter values. Module 30 can be 
programmed using Matlab™ software supplied by The Math Works Inc to implement 
network optimization methods. Optimization module 30 is linked to the experimental data 
entry module 16 to transfer thereto the generated set of optimal parameter values, which 
module 16 also stores the actual property values obtained from an experimental run 
based on the set of optimal parameter values. All experimental data is then transferred 
to the evaluation module 18 as mentioned before. 

A preferred embodiment of an optimization method according to the present 
invention will now be described with reference to Figs. 1 and 2. As illustrated in the 
general flow chart shown in Fig. 2, the method comprises a first step 40 of assigning 
values to a set of k property weights w ; representing relative importance of the k 
properties Y 3 for the characterization of the product, which properties are likely to be 
affected by the parameters of the process, from a modeling of the problem expressed as 
a hierarchical tree of these properties. Initial modeling and weight value generation are 
preferably performed using a method known as analytic hierarchy process (AHP), which 
was first proposed by t.W. Saaty, and more recently described in "Using The Analytic 
Hierarchy Process For Decision Making In Engineering Applications: Some Challenges" 
Triantaphyllou et al, International Journal of Industrial Engineering: Application and 
Practice, Vol. 2, No. 1. pp. 35-44, (1995), which is incorporated herein by reference. 

The AHP method consists of building a hierarchical tree from all properties, with 
one or more hierarchical levels depending on existing relations between the properties. 
For each level, a pair-wise comparison matrix is built between the properties of this level 
and presented at an input of the parameters weighting module 14 shown in Fig. 2, which 
executes in step 40. For each pair-wise comparison, the normalized eigenvector is 
derived associated with the higher eigen value. The components of this eigenvector give 
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the relative importance of each property called the local weight. Finally, the above 
normalized vectors are combined to find the global weight for each property. 

In a parallel direction, each pair-wise comparison is associated with a consistency 
index reflecting the transitivity relation between all comparison by pairs given by the 
formulator. Multi-criteria analysis software which is commercially available, such as 
Expertchoice™, Criterium™ or Ergo™, may be used to program module 14. For 
example, to one or more m main properties classified at a first (higher) level, may 
correspond one or more groups of properties classified at a second (lower) level, the 
latter properties being therefore identified as sub-properties. For each main property 
associated with a group of p sub-properties, a matrix of dimension Lp + lx/? + l] is built 
and filled, as a result of a pair-wise comparison between each property and sub- 
property, using relative importance values selected from a standard AHP scale. Next, a 
suitable algorithm performed by parameter weighting module 14 consists of first 
calculating the higher eigenvalue of the resulting numerical matrix, and then deriving a 
normalized relative importance vector of dimension [p + l] by an estimation of the left 
principal eigenvector of that matrix associated with the calculated main eigenvalue of the 
input matrix. The above algorithm is then applied to compare the m main properties of 
the higher level, from a pair-wise comparison matrix of dimension [wixm] from which a 
normalized relative importance vector of dimension [m] is derived. Finally, the above 
normalized vectors are combined according to the hierarchical relations to generate a 
global relative importance weight vector for the k properties of dimension + or 
[k]. In practice, it is generally appropriate to retain only each group of sub-properties 
without the corresponding main property, the sum of the weights related to the retained 
k properties/sub-properties being always equal to unity. 

According to the next step, namely step 42, parameter data and property data 
values are provided, which data is obtained from experimental runs using different sets 
of parameter values for the process, the various values for each parameter being 
preferably selected according to an expected operation range within which an optimal 
parameter value is likely to be found. The parameters X t used in the experimental runs 
should cover the extremes of the expected operational range for each parameter. 
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Generally, the number of formulation combinations required to determine an optima! 
formulation depends on many factors among which the more important ones are: 1) the 
formulation designer experience; 2) complexity of the formulation; 3) the availability of 
literature and experimental data available on the desired product; and 4) the analytical 
5 laboratory workload and throughput. According to the method of the present invention, 
the minimal number of experimental runs / to perform has been found to be equal to 
Ti + l , wherein n is the number of relevant parameters involved. A greater number of 
runs is certainly possible. Step 42 is performed by experimental data entry module 16 
shown in Fig.1. 

10 The method then comprises a step 44 of establishing property behavior mathematical 
relations linking the properties with the parameters and interactions thereof, in 
polynomial form. These property behavior relations provide an estimated property Ye } 
2* for each of the k properties Yj in terms of a number n of parameters X t from the 

parameter data and associated property data provided at step 42. Step 44 is typically 
I135 comprised of four sub-steps, namely 1) a parameters reduction step performed by 
jfi module 22, 2) a parameters interaction analysis step performed by module 20, 3) a 
^ property behavior modeling step performed by module 24, and 4) a property behavior 
Q relations generating step performed by module 26, as shown in Fig. 1. As to sub-step 1), 
% to provide a more efficient algorithm, from an initial number of identified parameters, the 
^0 most significant parameters, i.e. those significantly affecting each property, are identified 
to generate a reduced number n of significant parameters, especially where the initial 
number of identified parameters is greater than 8, as mentioned before. For that 
purpose, a statistical analysis algorithm can be used, which is based on parameter 
correlation calculations using parameter and property experimental data provided at 
25 prior step 42. Having obtained data related to / experimental runs involving an initial 
number p of parameters and a number k properties Y it each correlation factor 
contained in the correlation matrix is retained as significant whenever it complies with a 
predetermined condition in the following form: 

a < < b or - c < p g < - d (1 ) 

30 wherein a 9 b,c and d are predetermined limit values, typically set as follows: 
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0.5 < p & < 0.95 or -0.95 < p f < -0.5. (2) 

The parameters associated with the retained correlation factors form the reduced set of 
n parameters. 

It can be also shown that a minimum number / of runs at least equal to n + l is 
required to obtain reliable parameters estimation. Then, parameter interactions, that 
are in the form X ( Xj with i * j and which are significant, can be identified using the 
above relations (1), with the suggested specific ranges given in (2). The values for X t 
from the / experimental runs are combined with the retained correlation factors p & to 
form a final matrix W f with each element of the first column being equal to unity for the 
purpose of following sub-step 4). As to sub-step 3), it consists of establishing, for each 
property y. , a best model in terms of retained parameters and parameter interactions. A 

standard variance analysis is carried out to confirm relevancy of all parameter 
coefficients and parameter interaction coefficients, and to select by successive variance 
analysis operations through the use of modules 24, 20 and 22, a suitable model 
amongst different predetermined models of upgraded degrees, whenever difference in 
performance between a given model of degree r and a following model of degree r+1 is 
found to be not significant. The resulting best model is taken along with matrix Wand 
property experimental data in matrix y, as inputs for following sub-step 4) aimed at 
generating property behavior relations for each property Yj . A matrix C of coefficient 

values is given by the matrix: 

C = (w T w)' 1 xW T xY (8) 
having a dimension of [zw,*], wherein iw = n + / + l, t being the number of parameters 
interactions X i X j . Hence, estimated property values are given by: 
*i = fi(X 19 ... 9 X n >.„,X i Xj) 



Ye =C T X = < 



Y J =f j (X X9 ... 9 X n ^ 9 X l X j ) (9) 

Y k = f k (Xl,... 9 X n9 :., XfXj) 
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A following step 46 as shown in Fig. 2 aimed at generating a goal function is carried out 
by the module 28 shown in Fig. 1, from the set of k property weights w y produced at 
step 40, from the property behavior relations produced at step 44 and from the specified 
goal values for the properties Y r The basic goal vector can be expressed as follows: 

wherein O t is the specified goal values for the properties Y n with i = 1,...,* . 
The goal function to be minimized may be expressed as follows: 

(11) 

which goal function is expressed in terms of property weighted deviations 
between estimated values Ye } for the properties Y } and corresponding specified goal 
values 0 { for the same properties Y r A next step 48 as shown in Fig. 2 therefore 
consists of minimizing the goal function G to obtain a set of optimal parameter values 
for the parameters X if which step 48 is performed by module 30 shown in Fig. 1. 
Optimization step 48 generally can consider constraints on the parameter values in the 
form of one or more ranges, typically in a form (a i9 b t ) wherein a f <X g < b f , within which 
optimal parameter values shall be found, according to the type of each variable (i.e. 
discrete such as binary values, or continuous). 

The "G" goal function is determinated by experimentation. The optimization of the "G" 
function is a step by step procedure. The first step is to obtain the behavior laws with 
the best fit between the experimental data and their corresponding ideal value factor. 

The second step, the optimization is based on a initial point. 

x 0 = lx? 9 x 0 2 , ,x°„] 

and 

X k ~ l =X k -aHg\ x xk 
where 
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g=*G = gradient G 
H= the Hessian of G 

And based on the following goal function, we use the dimension reduction method by 
successive iterations 

1=1 

these iterations passed by 
/Or,.....*.-.) =0 

if i = \,....,k-\ 

and 

X n = f*(Xw-*Xn-l) 

Now the goal function can be: 

G{X x ,...,X n ) = G{X x ,...,X n _J- k {Xi ) 

We observe a perfect overlap between the two goal functions and on the stationary point 
the goal function will be; 



3G(X r ) _ 
dX, 



dG(x) dG(X)8f k (X lt ..JC n J 

ax. ex. ex, 



-I x K {x t t —,X„_ 1 ) 



These equations supply the maxima and minima of the goal functions including the 
maxima and minima from the starting goal function. 

This mathematical approach induces a reduction of the dimension of the variables, 
consequently we pass from "n" variables to "n-1" variables. In the actual case, we start 
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with the most important variables from the behavior laws with the highest weight values 
of the factor. 

This approach is known under the name of network optimization, in this case the 
5 network nodes are built by the optimal values of the variable by decreasing order of the 
factor's rank. 



After the iterative optimization step 48 is completed, although the set of optimal 
10 parameter values X0 { obtained can generally be considered as the solution to 
recommend, that solution is preferably evaluated amongst other alternative solutions by 
following steps 50 and 52 as shown in Fig. 2. At step 50, an experimental run of the 
2 process is carried out using the obtained set of optimal parameter values, to obtain 

I s * experimental values for the k properties Y x . The optimal property values X0 ( and 

Ul 

Fd 5 associated experimental property values are then evaluated at step 52 to obtain ranking 

*S 

jjj thereof amongst a number m of other alternative solutions, which may have been 

^ selected from knowledge base 12 shown in Fig. 1. This evaluation is preferably 

O performed by a complete AHP process algorithm, using the set of k property weights 

p Wj as previously obtained through step 40. 



Concrete Application 



An example illustrating an application of the method according to the present 
invention in the pharmaceutical field will now be described. 
25 Formulation and production process for enalapril maleate tablets were optimized 

in order to provide a drug product with satisfactory biological performance as well as 
stability when packaged and stored under ICH (International Conference on 
Harmonization) conditions. Three (3) independent formulation and process parameters 
(n = 3) were identified as having an impact on the stability of the drug product: 1) the 
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degree of drug neutralization during granulation (X^) ; 2) the manufacturing technology 
(X 2 ) ; and 3) The drug-to-excipient ratio in the formulation, i.e. dose strength (X^) . 

As to the degree of drug neutralization during granulation {X x ), it was classified 
as either complete, partial or no neutralization. In the case of complete neutralization, 
5 the drug and the alkaline agent were both added to the granulation fluid, i.e. water. 
Therefore, the alkaline agent neutralized the drug prior to its addition to the powder 
blend for the granulation procedure. In partial neutralization, both the drug and the 
alkaline agent were added to the powder blend, blended and water added as the 
granulation fluid for the granulation procedure. When water and/or alkaline agent were 
10 not added to the formulation, the drug was not neutralized. The level of water added as 
well as the drug-to-alkaline agent ratio were kept constant for all of the formulations. The 
m level of the alkaline agent was determined by the stoichiometry of the reaction. 
P The manufacturing technology {X 2 ) was either wet granulation (X 2 = 0 ) or direct 

jap 

!J! compression (X 2 =l). These two technologies are used worldwide for the 

W 

.pS manufacturing of probably more than 90% of all of the solid oral dosage forms. In the 

in 

(Tj wet granulation technology, the drug and other functional materials added to impart 
good processing attributes to the drug, often called excipients, are first blended together 
ifl and agglomerated into larger particles by the addition of a granulating fluid. The role of 
\f % the granulating fluid is to promote the development of adhesive forces between the 
*f£0 materials required for the agglomeration process. After granulation, the granulating fluid 
is removed by drying. When a direct compression approach is selected as a 
manufacturing method, the drug is first blended with the excipients and tablets produced 
without the use of a granulating fluid. 

As to the dose strength (X z ) t four doses of the product were developed, which 

25 were obtained by using two formulations with different drug-to-excipient ratios 
(continuous parameter values) compressed at different tablet weights. 

A total of nine (9) experimental runs involving different formulations based on a 
combination of the three parameters were prepared, as shown in Table 1. 



Run 
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1 


1.73 


1 


3.33 


2 


1.73 


0 


3.33 


3 


1.73 


o 


3.33 


4 


o 


o 


3.33 


5 


0 


1 


3.33 


6 


6.9 


1 


13.3 


7 


6.9 


0 


13.3 


8 


6.9 


0 


13.3 


9 


0 


0 


13.3 



TABLE 1 



p The nine formulations covered ail of the six (6) possible combinations for the wet 
granulation technology and three (3) combinations of direction compression. Tablets 
rU were manufactured by using enalapril maleate with USP/NF and EP excipients. In the 
$11 direct compression technology, there is not a sufficient amount of moisture to dissolve all 
v the drug and alkaline agent and provide for any significant neutralization reaction. 
O However, excipients do contain a certain level of adsorbed free moisture capable of 
g40 creating a microenvironment where small quantities of the drug and alkaline agent can 
^ be dissolved and become available for the neutralization reaction. These phenomena 
ifl could be responsible of the appearance of physical as well as chemical stability 
problems and where taken into account by evaluating three (3) formulation 
combinations. The nine (9) formulation combinations where prepared and the tablets 
15 were stored in opened containers at 25°C/60%RH and 40°C/75%RH for a 2-week 
period. These open container studies are typically conducted during the early 
formulation development phases of a product to purposely accelerate physical and 
chemical changes in formulations in order to select the lead candidate, i.e., the 
formulation with the best stability profile. After the 2-week time period, the tablets were 
20 removed from the environmental chambers and sent to the analytical department for 
their performance evaluation. The performance of the formulations was determined by 
measuring ten (£=10) properties as a function of time and temperature, which properties 
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were selected as follows, according to a hierarchical tree comprising properties and sub- 
properties: 

- Y11, Y12, Y13: % drug dissolved at 5, 15, and 30 min. (sub-properties of Y1); 

- Y2: % of cyclization product at time zero; 

- Y31, Y32: % cyclization product after 2 weeks at 25°C/60%RH and at 
40°C/75%RH (sub-properties of Y3); 

- Y4: differential between theoretical and actual assay in mg at time zero; 

- Y51 , Y52: differential between theoretical and actual assay in mg after 2 
weeks at 25°C/60%RH and at 40°C/75%RH (sub-properties of Y5); 

- Y6: % hydrolytic product after 2 weeks at 40°C/75%. 



Applying the AHP process with the standard scale for these properties, the decision 
matrixes given in Table 2 for the properties and in Tables 3, 4 and 5 for the sub- 
properties were built. 



Goal 


Y1 


Y2 


Y3 


Y4 


Y5 


Y6 


Y1 


1.0 


(5.0) 


(5.0) 


3.0 


1.0 


(5.0) 


Y2 




1.0 


(5.0) 


1.0 


(3.0) 


(7.0) 


Y3 






1.0 


7.0 


1.0 


1.0 


Y4 








1.0 


(7-0) 


(7.0) 


Y5 










1.0 


(1.0) 


Y6 












1.0 



TABLE 2 



Y1 


Y11 


Y12 


Y13 


Y11 


1.0 


7.0 


9.0 


Y12 




1.0 


3.0 


Y13 






1.0 
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Y3 


Y31 


Y32 


Y31 


1.0 


(5.0) 


Y32 




1.o 


TABLE 4 


Y5 


Y51 


Y52 


Y51 


1.0 


(5.0) 


Y52 




1.0 



TABLE 5 



From the decision matrixes, the following weight values for the £=10 properties/sub- 
properties are given in Table 6, the sum of the weights being equal to unity. 



(Y11) 


w 2 
(Y12) 


w 3 
(Y13) 


^4 

(Y2) 


(Y31) 


w 6 
(Y32) 


W 7 

(Y4) 


(Y51) 


w a 
(Y52) 


W 10 

(Y6) 


0.060 


0.011 


0.005 


0.098 


0.046 


0.230 


0.039 


0.034 


0.171 


0.305 



TABLE 6 



Experimental property data that were obtained from nine (9) runs of the process using 
the selected nine (9) combinations of parameter values of Table 1 , are given in Table 7. 



Run 


Y11 


Y12 


Y13 


Y2 


Y31 


Y32 


Y4 


Y51 


Y52 


Y6 
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1 


90.40 


102.7 


102.1 


0.88 


1.210 


15.06 


0.080 


(0.01) 


0.660 


4.260 


2 


91.20 


95.10 


95.80 


0.63 


1.070 


3.510 


0.120 


0.020 


0.340 


5.350 


3 


86.80 


100.4 


101.6 


0.90 


0.850 


0.720 


0.080 


0.030 


0.110 


3.670 


4 


94.00 


96.70 


97.20 


0.85 


1.340 


16.88 


0.300 


0.160 


0.780 


2.590 


5 


75.40 


106.2 


107.1 


0.83 


1.350 


22.54 


(0.08) 


0.040 


0,960 


3.570 


6 


94.70 


98.20 


98.30 


0.54 


1.110 


8.600 


0.380 


0.270 


2.200 


6.730 


7 


90.20 


98.70 


98.40 


0.64 


1 .670 


b.you 


U.O IU 




1 7^0 


7.710 


8 


51.70 


100.4 


100.8 


0.50 


0.75 


0.750 


(0.04) 


(0.14) 


0.310 


3.580 


9 


84.90 


95.30 


95.20 


0.55 


0.11 


0.110 


0.930 


0.740 


2.210 


1.450 



TABLE 7 

Since n = 3 < 8, the parameter reduction step is not required for the purpose of the 
instant case. As to the statistical analysis of parameters interaction, since a correlation 
factor /7 13 =0.7013 for the *,* 3 interaction was calculated, that interaction can be 
considered as significant since the condition 0.5 < /? 13 < 0.95 is satisfied. The following 
property behavior relations were established: 



Y11 =81.916+4.56^ + 4.074 * 2 + 0.224 * 3 -0.423*,* 3 ; 
Y12 = 101.93 -1.45 X, + 3.81 * 2 -0.51* 3 +0.14X l X 3 ; 
Y13 = 102.16 -14.8*, + 3.46* 2 -0.52* 3 + 0.14*,* 3 ; 
Y2 = 0.92 - 0.025 *, + 0.025 X 2 - 0.028 X, + 0.0018 *,* 3 ; 
Y31 = 1.42-0.23*, + 0.057 X 2 - 0.03*3 + 0.019*,*3; 
Y32 = 17.18 - 8.78*, + 8.15* 2 - 0.46X 3 + Q.56*^; 
Y4 = -0.15 + 0.0193*, - 0.022 * 2 + 0.08X, - 0.0091 *,* 3 ; 
Y51 = -0.135 - 0.028*, + 0.0385* 2 + 0.066* 3 - 0.0045*,* 3 ; 
Y52 = 0.00089 - 0.256*, + 0.63 X 2 + 0.166* 3 + 0.008^*3 ; 
Y6 =3.24 + 0.9*. +0.57* 2 - 0.13* 3 - 0.02 *,*,; 



(18) 
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The specified goal values for the properties as given in Table 8 were used to establi 
the goal function that was minimized to generate the following set of optimal parametei 



X, = 3.39 

X 2 = 0 (wet granulation) 
X 3 = 7.46 

The associated experimental property values are given in Table 9. 





o 2 


0 3 


o 4 


o 5 


o 6 




o % 


o, 


o 10 


(Y11) 


(Y12) 


(Y13) 


(Y2) 


(Y31) 


(Y32) 


(Y4) 


(Y51) 


(Y52) 


(Y6) 


100.00 


100.00 


85.00 


0.500 


0.500 


0.500 


0.050 


0.050 


0.050 


1.00 






TABLE 8 






Y11 


Y12 


Y13 


Y2 


Y31 


Y32 


Y4 


Y51 


Y52 


Y6 


97 


104 


106 


0.49 


0.58 


0.73 


0.031 


0.008 


0.047 


<0.05% 



TABLE 9 

Applying the method for the particular case where only the minimum four 
(n + 1 = 3 + 1 = 4) experimental runs required were used, runs 1, 3, 6 and 9 were 
selected to provide the parameter and property data as given in Table 7. As to the 
statistical analysis of parameters interaction, since a correlation factor p l3 =0.332 for the 
X x X y interaction was calculated, that interaction cannot be considered as significant 
since 0.5 < p l3 < 0.95 is not satisfied. The following property behavior relations were 
established: 
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Y11 =85.36 + 0.99*, + 3.60* 2 -0.34*,; 
Y12 = 101.90 + 0.086*, + 2.30 * 2 - 0.50 *,; 
Y13 = 102.87 + 0.38*, + 0.506 * 2 - 0.58*,; 
Y2 =1.013 + 0.001*, -0.019* 2 - 0.034 * 3 ; 

Y31 = 0.88 - 0.038*, + 0.36 * 2 + 0.0095*, ; (20) 

Y32 = 2.88 - 2.425*, + 14.33 * 2 + 0.61 *,; 

Y4 = -0.02 - 0.080*, - 0.0007 * 2 + 0.071 *, ; 

Y51 = -0. 0568 - 0.065*, - 0.02 * 2 + 0.06*,; 

Y52 = -0.4 - 0.08*, + 0.55 * 2 + 0.2 * 3 ; 

Y6 =2.84 + 0.68*, +0.59* 2 - 0.1*,. 

The same specified goal values for the properties as given in Table 8 were used to 
establish the goal function that was minimized to generate the following set of optimal 
parameters: 

*, = 3.32 

* 2 = 0 (wet granulation) (21 ) 

*, = 7.09 

The associated experimental property values are given in Table 10. 



Y11 


Y12 


Y13 


Y2 


Y31 


Y32 


Y4 


Y51 


Y52 


Y6 


78 


103 


106 


0.51 


0.55 


0.73 


0.034 


0.04 


0.029 


<0.05% 



TABLE 10 

Comparing the set of parameter values given at (20) with the former set obtained from 
all nine (9) experimental runs given at (19), it can be noted that both sets are very 
similar. Actually, from a pharmaceutical standpoint, they could almost be considered as 
identical. 
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What is claimed is: 

1. A method of optimizing parameter values in a process of producing a product, 
said process being essentially controlled by a set of n parameters X t . affecting a 

set of k properties Y. characterizing the product, said method comprising the 

steps of: 

i) assigning values to a set of k property weights w. representing 
relative importance of said properties T} for the characterization of said 
product; 

ii) establishing property behavior mathematical relations giving an 
estimated property Ye } for each said property Y } in terms of said parameters 

X t from given parameter data and associated property data; 

iii) using said property weights Wj to establish a goal function in terms 
of property weighted deviations between the estimated properties Yej and 
corresponding specified goal values for said properties 7 y ; and 

iv) minimizing the goal function to generate a set of n optimal 
parameter values for said parameters X t . 

2. A method according to claim 1, wherein said product is a composition, said set of 
optimal parameter values characterizing an optimal formulation for the 
composition. 

3. A method according to claim 1, wherein said product is a pharmaceutical product, 
said set of optimal parameter values characterizing an optimal formulation for the 
pharmaceutical product. 

4. A method according to claim 1, wherein the values for said property weights Wj 
are obtained using an algorithm based on an analytic hierarchy process. 
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5. A method according to claim 4, wherein said given property data are obtained 
through a number / of experimental runs of said process using said given 
parameter data, each said run using a distinct set of values for said given 
parameter data. 

6. A method according to claim 5, wherein / > n + 1 . 

7. A method according to claim 1, wherein said goal function is expressed as 
follows: 

8. wherein 0 } are said specified goal values for said properties Y r A method 
according to claim 7, wherein said minimizing step is performed by successive 
iterations of: 

G(^,...,^J = XU(^,..,^)r 

9. A method according to claim 8, wherein said goal function is minimized according 
to one or more specified ranges (a t ,b g ) wherein a t <X t .< b t for one or more of 
said optimal parameter values. 

10. A method according to claim 1 , further comprising the steps of: 

i) performing experimentally said process using said set of optimal 
parameters values to obtain corresponding experimental values for said 
properties 7/, 

ii) ranking said set of optimal parameters values over predetermined 
alternative sets of parameters values for said X % . 
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11. A method according to claim 10, wherein said ranking step is performed using an 
algorithm based on an analytic hierarchy process. 

12. A method according to claim 10, further including the step of: 

i) incorporating said set of optimal parameters values and said 
corresponding experimental values for said properties Yj respectively into 

said given parameter and associated property data; 

ii) repeating said steps ii) to iv) to generate a new set of optimal 
parameters values for said parameters X. . 

13. A method of producing a pharmaceutical product using optimized process 
parameter values, said process being essentially controlled by a set of n parameters X t 
characterizing a formulation for said product, said parameters X s affecting a set of k 
properties Y } characterizing the product, said method comprising the steps of: 

a) conducting a number of / of experimental runs of said process each using 
a selected distinct set of values for said parameters X i covering substantially all 
extreme values within a chosen range of values for each one of said parameters 
X i , wherein / is at least equal to n + 1 and is substantially less than a number 
used in the Fractional Factorial Matrix method; 

b) measuring values for said properties Y } characterizing the product in each 

of said / experimental runs, whereby parameter data and associated property 
data are obtained from said selected distinct set of values for said parameters X t 

and said measured values for said properties Y } , respectively; 

c) determining an importance of said properties Y } for the characterization of 
said product, comparing said importance of said properties Y } relative to one 
another, and assigning values to a set of k property weights Wj representing a 
relative importance of said properties Y. for the characterization of said product; 
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d) calculating a set of optimal parameter values for said parameters X i using 
said measured values for said properties Y } and said assigned values of said set 
of k property weights Wj ; and 

e) producing said pharmaceutical product using said optimized process 
parameter values X s calculated in the previous step. 

14. a method according to claim 13, wherein said step of calculating 

comprises: 

i) establishing property behavior mathematical relations giving an 

estimated property Ye } for each said property Yj in terms of said 
parameters X i from said parameter data and associated property data; 

jj) using said property weights w y to establish a process goal function 

in terms of property weighted deviations between the estimated properties 
Ye j and corresponding specified goal values for said properties I}; and 

iii) minimizing the process goal function to generate a set of optimal 

parameter values for said parameters X i . 

15. A method according to claim 14, wherein the values for said property weights w, 
are obtained by an algorithm based on an analytic hierarchy process. 

16. A method according to claim 14, wherein / = n + 1. 

17. A method according to claim 14, wherein said goal function is expressed as 
follows: 

GiX^.XJ^MjiYe^Ojf 
wherein O y are said specified goal values for said properties 7, . 
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18. A method according to claim 17, wherein said minimizing step is performed 
through successive iterations. 

19. A method according to claim 18, wherein said goal function is minimized 
according to one or more specified ranges (a„&,) wherein a t <X t < b, for one or more of 
said optimal parameters values. 

20. A method according to claim 14, further comprising the steps of: 
performing experimentally said process using said set of optimal parameters 
values to obtain corresponding experimental values for said properties Yy, 
ranking said set of optimal parameters values over predetermined alternative sets 
of parameters values for said X t . 

21. A method according to claim 20, wherein said ranking step is performed 
through an algorithm based on an analytic hierarchy process. 

22. A method according to claim 21 , further including the steps of: 

incorporating said set of optimal parameters values and said corresponding 
experimental values for said properties Y } respectively into said given parameter 

and associated property data; 

repeating said steps a), b) and d) to generate a new set of optimal parameters 
values for said parameters X t . 
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ABSTRACT OF THE DISCLOSURE 

A method of optimizing parameter values in a process for producing a product, which 
process is essentially controlled by a set of parameters affecting a set of properties 
5 characterizing the product. The method uses an analytic hierarchy process (AHP) to 
associate a weight with each property according to its relative importance to obtain 
desired product characteristics. The method also uses parameter data and measured 
property data from a required number of experimental runs of the process, from which 
data property behavior relations between each property and the parameters are 
10 statistically established, which relations give estimated property values. Using the 
property weights, a process goal function is established, which is expressed in terms of 
weighted deviations between the estimated property values and the corresponding goal 
values for the properties. Finally, the process goal function is minimized in order to 
2 generate a set of optimal parameter values for the process. The AHP process may be 
H5 further used to evaluate the optimal solution compared to a set of other alternative 
jij solutions. The method allows the use of input data obtained from a minimum number of 
jfj experimental runs, to generate a reliable set of parameter values. An example 
M illustrating an application of the method for the design of a pharmaceutical formulation is 

V. 

□ described. 
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Figure No. 2: 
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